arXiv: 1508.06381 vl [cs.IT] 26 Aug 2015 


1 


Joint Transceiver Design Algorithms for Multiuser 
MISO Relay Systems with Energy Harvesting 

Yunlong Cai, Ming-Min Zhao, Qingjiang Shi, Benoit Champagne and Min-Jian Zhao 


Abstract —In this paper, we investigate a multiuser multiple- 
input single-output (MISO) relay system with simultaneous wire¬ 
less information and power transfer (SWIFT), where the received 
signal is divided into two parts for information decoding (ID) and 
energy harvesting (EH), respectively. Assuming that both base 
station (BS) and relay station (RS) are eqnlpped with mnltiple 
antennas, this work studies the joint transceiver design problem 
for the BS beamforming vectors, the RS ampllfy-and-forward 
(AF) transformation matrix and the power splitting (PS) ratios 
at the single-antenna receivers. The aim is to minimize the total 
transmission power of the BS and the RS under both signal-to- 
interference-plus-noise ratio (SINR) and EH constraints. Firstly, 
an iterative algorithm based on alternating optimization (AO) and 
with guaranteed convergence is proposed to successively optimize 
the transceiver coefficients. This AO-based approach is then 
extended into a robust transceiver design against quantization 
errors in channel state information (CSI), by using semideflnite 
relaxation (SDR) and the S-procedure. Secondly, a novel design 
scheme based on switched relaying (SR) is proposed that can sig¬ 
nificantly reduce the computational complexity and overhead of 
the AO-based designs while maintaining a similar performance. 
In the proposed SR scheme, the RS is equipped with a codebook 
of permutation matrices. For each permutation matrix, a latent 
transceiver is designed which consists of BS beamforming vectors, 
optimally scaled RS permntatlon matrix and receiver PS ratios. 
For the given CSI, the optimal transceiver with the lowest total 
power consumption is selected for transmission. We propose 
a concave-convex procedure (CCCP) based and subgradient- 
type iterative algorithms for the non-robust and robust latent 
transceiver designs. Simulation resnlts are presented to validate 
the effectiveness of all the proposed algorithms. 

Index Terms —Beamforming, concave-convex procednre, 

switched relaying, power splitting, energy harvesting, SWIFT. 

I. Introduction 

Recently, electromagnetic (EM) energy transfer techniques 
have attracted considerable interest in the wireless research 
community. Indeed, by exploiting the radiative far-held prop¬ 
erties of EM waves, these techniques could in theory enable a 
radio receiver to harvest energy from its environment, thereby 
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relaxing the battery requirements on user devices. An interest¬ 
ing application of wireless energy transfer is to jointly transmit 
information and energy using the same waveform, which 
is known as simultaneous wireless information and power 
transfer (SWIFT). The idea of SWIFT was first proposed 
by Varshney in ||Tj, which characterizes the rate-energy (R- 
E) tradeoff in a discrete memoryless channel. The study of 
the R-E tradeoff was later extended to frequency selective 
fading channels in Q. In Q, the authors studied a two-way 
communication scenario with noiseless channels and limited 
resources, with emphasis on the tradeoffs due to the need 
of balancing the information flow with the resulting energy 
exchange among the communicating nodes. The authors of Q 
focused on multiuser systems and demonstrated that energy 
transfer constraints call for additional coordination among 
distributed nodes of a wireless network. However, the above 
studies have not been realized yet due to practical circuit 
limitations. 

The first practical receiver structure that makes SWIFT 
possible was proposed in Q, where two practical signal sep¬ 
aration schemes were considered, namely, the time switching 
(TS) scheme, where the receiver switches between information 
decoding (ID) and energy harvesting (EH), and the power 
splitting (FS) scheme, where the received signal is split into 
two streams, such that a fraction p (0 <P< 1) of the received 
signal power is used for ID while the remaining fraction (1—p) 
is used for EH. FS-based algorithms were considered in Q- 
|10|. In 1^, the authors derived the optimal FS rule at the 
receiver, for both SISO (single-input single-output) and SIMO 
(single-input multiple-output) systems, in order to optimize the 
R-E performance tradeoff. In Q, the authors proposed two 
practical receiver architectures, i.e., separated and integrated 
information and energy receivers, based on dynamic power 
splitting. These authors characterized the R-E performance by 
taking circuit power consumption into account. 

The authors of Q studied the joint beamforming and power 
splitting (JBFS) design for a multiuser multiple-input single¬ 
output (MISO) broadcast system with SWIFT. In this study, 
the total transmission power at the base station (BS) is mini¬ 
mized subject to signal-to-interference-plus-noise ratio (SINR) 
and EH constraints for all the receivers. The JBFS problem 
for a iT-user MISO interference channel was considered in 
0, where the authors used the semideflnite relaxation (SDR) 
technique to address the non-convex problem and proved that 
the SDR is tight in the case of AT = 2 or 3. Different 
from the approach of Q, an alternative second-order-cone- 
programming (SOCF) relaxation method was proposed in pO) . 
The SOCF relaxation based method guarantees a feasible so- 



2 


lution to the JBPS problem and has lower complexity than the 
SDR method, while achieving a performance extremely close 
to the minimum transmission power. A primal-decomposition 
based decentralized algorithm was also presented in this work. 

Besides, the SWIFT technique for relay systems was con¬ 
sidered in GD-EI. Specifically, the authors of GD proposed 
a joint source and relay precoding design algorithm to achieve 
different tradeoffs between the energy transfer and the infor¬ 
mation rate. In the relay beamforming design problem 
for the SWIFT scheme was considered in a non-regenerative 
two-way multi-antenna relay network, where a global opti¬ 
mal solution, a local optimal solution and a low-complexity 
suboptimal solution were proposed. The design of a high-rate 
beamformer that supports multiple communication pairs and 
intended for SWIFT in wireless relay networks was considered 
in p3| . In | |T4) , a game-theoretical framework was developed 
to address the distributed power splitting problem for SWIFT 
in relay interference channels. The aforementioned works 
assume the availability of perfect channel state information 
(CSI), which may be impractical in realistic implementations 
of these systems. 

In this paper, we consider the use of FS-based receivers in 
a general multiuser MISO relay systems, i.e., where a multi¬ 
antenna relay station (RS) is incorporated into the traditional 
MISO multiuser setup, as in GD-lIZ)- We focus our study 
on the downlink transmission where the BS first transmits 
the signals intended for different receivers via beamforming 
to the RS. Then, the RS processes the received signals 
through an amplify-and-forward (AF) transformation matrix 
and broadcasts it to all the receivers. Finally, each receiver 
employs the FS technique to decode information and harvest 
energy simultaneously. We present an optimization framework 
for the joint design of the BS beamforming vectors, the RS 
AF transformation matrix and the receiver FS ratios aiming to 
minimize the total power consumption under a set of minimum 
SINR and EFI constraints at the receivers. In this work, we 
shall assume that the proposed joint design algorithms are 
implemented at the BSQ 

To this end, we first propose an iterative algorithm based 
on alternating optimization (AO) and with guaranteed con¬ 
vergence to successively optimize the transceiver coefficients, 
i.e. the BS beamformers, the AF transformation matrix and 
the receiver FS ratios. We show that each subproblem can 
be relaxed as a semidefinite programming (SDF) problem by 
applying the celebrated SDR technique m- This AO-based 
approach is then extended into a robust transceiver design 
against quantization errors in CSI based on the SDR technique 
and the S-procedure | |22| , p^ . While the performance of 
the AO-based designs is remarkable, we note that in this 
approach, the BS needs to compute and send the complete 
AF matrix to the RS before transmission, which entails 
high design complexity and signaling requirements for the 
overhead. Secondly, to reduce the design and implementation 
costs, we present extensions of the switched relaying (SR) 

’in practice, in cellular systems it is preferable to implement most of the 
signal processing operations at the BS rather than the RS due to the fact that 
the BS is more powerful and the RS is expected to have a simple structure 
and low energy consumption |18|-|20|. 


processing reported in p4) to multiuser MISO relay systems 
with energy harvesting. For this scheme, we equip the RS 
with a codebook of permutation matrices. Based on each 
element of the codebook, we can obtain a permuted channel 
matrix and therefore create a latent transceiver, which includes 
the BS beamforming vectors, an optimally scaled relaying 
permutation matrix and the receiver FS ratios. Among the 
latent transceivers so obtained for given CSI, the optimal one 
is chosen according to a suitable criterion for transmission. 

Specifically, SR-based algorithm constructs the RS AF 
transformation matrix within each latent transceiver by mul¬ 
tiplying a permutation matrix from the codebook with a 
power scaling factor. Before data transmission, the BS sends 
the index of the permutation matrix (instead of the RS AF 
transformation matrix), the RS power scaling factor, and the 
FS ratios corresponding to the optimal transceiver to the 
RS and the receivers through signaling channels. Compared 
to the AO-based algorithms, the SR approach reduces the 
number of optimization variables from the complete RS AF 
transformation matrix to a single power scaling factor, which 
in turn significantly reduces the computational complexity 
and signaling overhead. To design the latent transceivers in 
the SR scheme, an iterative algorithm based on the concave- 
convex procedure (CCCF) p5) , |26| is proposed that uses 
the estimated CSI, and is guaranteed to converge to a local 
optimal point. Each subproblem in the iterative algorithm can 
be formulated as an SOCF problem. By taking the quantization 
CSI errors into account, we propose a robust subgradient algo¬ 
rithm which utilizes the side information provided by standard 
convex solvers. Eurthermore, a simplified SR-based transceiver 
design algorithm is proposed with much less computational 
complexity, while achieving a performance extremely close 
to the non-simplified scheme. Einally, two efficient codebook 
design approaches are developed. 

Simulation results demonstrate that the proposed transceiver 
design algorithms are capable of providing robustness against 
the effects of norm-bounded CSI errors. In particular, the 
SR-based algorithms using a small codebook of permutation 
matrices can achieve almost the same performance as the 
AO-based algorithms but with reduced design/implementation 
complexity and signaling overhead. 

The reminder of this paper is organized as follows. Section 
[n| presents the multiuser MISO relay system model, the 
channel error model and the problem formulation. In Section 
[nl| the AO-based transceiver design algorithms for both non- 
robust and robust cases are developed. Section IV discusses 
the SR-based transceiver scheme, including the proposed la¬ 
tent transceiver design algorithms and the codebook design 
methods. A complexity analysis of the proposed algorithms 
along with a discussion of their initialization are provided in 
Section [V] Einally, in Section computer simulations are 
used to verify the proposed algorithms. Conclusions are drawn 
in Section l\^ 

Notations: Scalars, vectors and matrices are respectively 
denoted by lower case, boldface lower case and boldface upper 
case letters. Eor a square matrix A, Tr(A), rank(A), A^, 
conj(A), and A^ denote its trace, rank, transpose, conjugate, 
and conjugate transpose respectively, while A A 0 means that 
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A is a positive semi definite matrix. The operator vec(-) stacks 
the elements of a matrix in one long column vector, invp(a;) 
denotes the inverse of the positive portion, i.e., . |j • ||, 

(•)!, and | • | denote the Euclidean norm of a complex vector, the 
factorial operator, and the absolute value of a complex scalar, 
respectively. denotes the space of m x n 

complex (real) matrices, and K+ (M_) denotes the set of 
positive (negative) real numbers. Finally, the symbol ® denotes 
the Kronecker product of two vectors/matrices. 


11. System Model and Problem Formulation 
A. Proposed System Model 




BS 




w 





Fig. 1. The multiuser MISO relay system with PS-based receiver. Each 
receiver splits the received signal into two parts for ID and EH, respectively. 


In this work, we consider a multiuser MISO relay system 
which consists of one BS, one RS, and K mobile receivers 
indexed byfcs/C = {l,...,iT}. The number of antennas at 
the BS and the RS is denoted as Nt and Nj., respectively, while 
each receiver is equipped with a single antenna. We assume 
that K < mm{Nt, Nr}, which provides sufficient degrees of 
freedom for signal detection. We also assume that the classical 
two-hop AF relaying protocol is employed, and that the 
direct links between the BS and the receivers are sufficiently 
weak (to be ignored). Different from the conventional mul¬ 
tiuser MISO relay channels | [T7| , we here consider PS-based 
receivers. The received signal at each receiver is split into 
two separate signal streams with different power levels; one is 
sent to the EH receiver and the other one is diverted to the ID 
receiver Q, as shown in Fig.[T] Under the above assumptions, 
the signals are transmitted in two phases as explained below. 

In the first phase, the BS transmits K data streams, each 
carrying an independent message intended for one of the K 
receivers. Thus, the transmitted data vector at the BS can be 
expressed as 

K 

= ( 1 ) 

where Sk is the data signal for receiver k, with zero mean and 
variance Fjlsfel^} = 1, and fj. € denotes the transmit 

beamforming vector. The data signals assumed 

to be independent of each other. The transmit power of the 


BS can be shown as 

K 

PB = i?{xBxf} = 5]||ffef. (2) 

k=l 

The transmission from the BS to the RS can be modeled as 
a standard point-to-point MIMO channel. Hence, the received 
data vector at the RS can be expressed as 

K 

yi? = + n^, (3) 

k^l 

where G S £NrxNt denotes the MIMO channel from the BS 
to the RS, and the complex circular Gaus¬ 

sian noise vector at the RS, with zero mean and covariance 
E{'n.r'n.^} = CT^I, where is the average noise power. It is 
assumed that the transmitted signals are independent 

of the noise vector Ur- 

In the second phase, the RS forwards the received signal to 
all the receivers after performing linear AF processing. Hence 
the vector signal transmitted from the RS can be formulated 
as 

Xfl = WyR, (4) 

where W S C^r-xNr transformation matrix at the 

RS. The transmission power of the RS can be shown as 

K 

Pr = Fjxfixf} = ^ IIWGffef + a2||W||2. (5) 

k^l 

Finally, The signal received at the kih receiver. A: € /C, is given 
by 

K 

2/fe = hfWGy^ffcSfc-f hfWrir+ nfe, (6) 

fc=i 

where hj. S denotes the complex conjugate channel 

vector between the RS and receiver k, and is the additive 
noise introduced by the receive antenna at receiver k, which 
is assumed to be a complex circular Gaussian variable with 
zero mean and variance a^. 

Fet Pfe (0 < pfc < 1) denotes the PS ratio for receiver k, 
which means that portion pk of the signal power is used for 
signal detection while the remaining portion 1 — Pk is diverted 
to an energy harvester. Thus, on the one hand, the signal 
available for ID at receiver k can be expressed as 

r® = WG ikSk + hf Wn, + + Vk, (7) 

where Vk is the additional complex circular Gaussian circuit 
noise with zero mean and variance E{\vk\'^} = oJj., resulting 
from phase offsets and non-linearities during baseband conver¬ 
sion Q. Thus, the SINR at the fcth receiver can be expressed 
as 

^ PfcIhfWGffcp 

i fc = -V-v- 

Pk [ £ |h|^WGf,|2 + a2||hfW|P + a2 ) +^2 

\jAk J 


( 8 ) 
























4 


On the other hand, the total harvested energy that can be 
stored by receiver k is equal to 

^r=6(l-Pfc) + a2||hfwf+ ^2^ , 

(9) 

where G (0,1] denotes the energy conversion efficiency of 
the /cth EH unit, which indicates that only a portion of the 
radio frequency energy can be stored. 

B. Channel Error Model 

We assume that the required CSI can be estimated at the 
BS and the RS by means of a suitable channel estimation 
algorithms Il2§-||^. In this work, we consider a time division 
duplex (TDD) mode and assume that the radio channels vary 
sufficiently slowly over time, so that the downlink transmit 
CSI can be obtained by channel reciprocity m- Therefore, 
the BS can obtain the MIMO channel matrix G involved the 
hrst phase of transmission, while the RS can obtain the channel 
vectors {h^} involved in the second. However, the second 
phase CSI needs to be fed back from the RS to the BS using 
signaling channels ll3T)-|[^ for implementation of the design 
algorithms, as explained earlier. Hence, here we assume that 
the BS can perfectly know the hrst phase CSI, matrix G, 
but not the second phase CSI, i.e. vectors {h^}. We consider 
quantization errors during channel feedback from the RS to the 
BS and employ the NBE model p4| to characterize this type 
of imperfection. In particular, the true (but unknown) second 
phase CSI can be expressed as follows: 

hfc = hfe + Gfe, k G IC, (10) 

where denotes the estimated channel vector, while 
denotes the CSI error vector. We assume that vector is 
bounded in its Euclidean norm, that is 

Ijefcll < 77fc, kG /C, (11) 

where is a known positive constant. Equivalently, 
belongs to the uncertainty set 3?^ dehned as 

= {h|h = hfc + Bk, ||efc|| < T]k}, kGiC. (12) 

The shape and the size of 3?fc model the type of uncertainty 
in the estimated CSI, which is linked to the physical phe¬ 
nomenon producing the CSI errors. It should be emphasized 
that the actual errors e/^ are assumed to be unknown while 
the corresponding upper bounds r]k can be obtained using 
preliminary knowledge about the type of imperfection and/or 
coarse knowledge of the channel type and main characteristics 



C. Problem Formulation 

In this subsection, we formulate the optimization problem 
for the joint design of W, {f/dfeex: {pk}kGK so as to min¬ 

imize the total power consumption at the BS and the RS under 
the constraint that a set of minimum SINR and EH targets be 


satished at the receivers. In this study, we consider both non- 
robust and robust designs (against quantization errors)]^ The 
non-robust optimization problem can be expressed a^ 

min Pb + Pr 

W, {ffc. pk} 

s.t. Tfe > 7 fc, P™ > (13) 

0 < Pfc < 1, Vfc e /c, 

where in the evaluation of (|^ and Q, the true channel vectors 
hfc are replaced by their estimates h^. 

Similarly, the robust optimization problem can be formu¬ 
lated as 

min Pb + Pr 

W. {ft, pk} 

s.t. Tfe > 7 fc, P™ > (14) 

0 < Pfc < 1, llefclp < T]l, yk G K., 

where in this case, the channel vectors in (|^ and (|^ are given 
by ([g. 

It is not difficult to see that both CD and d are in 
general non-convex because both their objective functions 
and constraints are not convex over W, {fj,} and {pk}- 
Eurthermore, d involves an inhnite number of constraints. 
These characteristics make it intractable to obtain the global 
optimal solution for ( |T3] l and ( [T4) i. In the sequel, we present 
two algorithms for obtaining suboptimal solutions to the above 
problems by applying proper convex optimization techniques, 
which are based on AO and SR, respectively. We complete 
this section with the following lemma. 

Lemma 1: Problem ( [T3] l is feasible for any hnite user SINR 
targets if rank(HG) = K, where H = [hi,..., h.^]^. 

Proof: The feasibility of d is not connected to the EH 
constraints and PS ratios according to Lemma 3.1 & Lemma 
3.2 in p^ . Hence, according to Theorem 1 in p7| , we can 
easily see that the above lemma holds. 

On the other hand, the feasibility of problem d has 
not been well studied in the literature and still remains an 
open question, which would be an interesting topic for future 
research. 

HI. Alternating Optimization Based Joint 
Transceiver Design 

In this section, we present the AO-based transceiver design 
algorithms for the jointly optimization of the BS beamforming 
vectors, the RS AE transformation matrix and the receiver PS 
ratios, for both non-robust and robust cases. In the proposed 
design, {f^, pk}k^K W are successively optimized in turn 
with the other hxed. We show that each subproblem for the 
optimization of {f^, pk} or W can be reformulated as an SDP 
problem based on the celebrated SDR technique and the S- 
procedure. Eurthermore, modihed randomization techniques 
based on the worst-case concept are provided to recover a 
rank-one solution when higher-rank solutions are returned. 

^The non-robust joint design is only based on the estimated CSI, while the 
robust design algorithm takes the CSI error into account. 

^There is in general a power consumption tradeoff between the BS and 
the RS dz). and the objective function should be Pq + olPr where a is 
a positive weight. We assume that a = 1 in this work, since this does not 
change the nature of the problem. 
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A. Non-Robust AO-Based Joint Transceiver Design 

In the following, we introduce a non-robust joint transceiver 
design algorithm. First, let us consider the optimization of 
{ffc,Pfc} in while the RS AF matrix W is fixed. The 
celebrated SDR technique can be applied to solve the re¬ 
maining optimization problem by introducing a new variable 
Ffe = ffcf^. Hence, problem can be reformulated as the 
following problem by ignoring the rank-one constraints for all 
Ffe's: 


min X; Tr(Ffc) 

{Ffc, pk} k=l 


E Tr(Q,) 


K 


S-t- - E hf Qj-hfc > 


K 


tlik 

ikii — Pk 


(Tk <J„ 


E hf Qjhfc > 

Ffc A 0, 0 < pfc < 1, Vfc G /C, 


IIWII 


?l|hfW|l2+a^ 


h?W| 


Pk ■ 


(15) 


where = WGFfcG^W^. In this way, problem (13i is 
relaxed to a convex SDP problem, which can be efficiently 
solved by off-the-shelf algorithms p7) . Let {F^} and {p^} 
denote the optimal solution to Based on Proposition 4.1 
in p^ , it can be verified that {F^} and {pj} satisfy the 
first two sets of constraints of problem ( [T5] l with equality and 
{F^} satisfy rank(F^) = 1, Vfc. Thus the optimal solution of 
problem ( [T3| ) can be expressed as {fj^jp^}, where is the 
principal component of F^, such that F^ = llffell = 

yjfk and fk is the largest eigenvalue of F^. 

Next, we consider the optimization of the RS AF matrix W 
while assuming that {fk,Pk} are fixed. Noting that x^Yz = 
vec(xz^)^vec(Y), the numerator in ^ can be rewritten as 


Pfc|hfWGffe|2 = p,|g^,vec(W)|2 

= PfcgfefcWconj(gfcfc), 


(16) 


where g^j = vec(conj(hfe)fJ^G^) and W = 
vec(W)vec(W)^. Similarly, the denominator in ^ 
can be reformulated as 

K _ 

Pk E grjWconj(gfcj) 

^ ^ (17) 

+Pk<Jr E hf EjWE^hfc -f pkal -f 

J=i 

where G {0, ^ linear mapping matrix such 

that h|^EjVec(W) = h^W(:,j), where W(:,j) denotes the 
jth column of W. With the help of ( [T6| and ( [T7| ), the SINR 
and EH constraints in problem ( [T3| ) can be expressed as the 
following two inequalities 


K 


Ik 


gLWconj(gfcfc) - X; gLWconj(gfcj)- 


j^k 


N-r 


(18) 


^ hfE^WEfh, > ^2 + 


K 


E gfe, Wconj(gfcj)-f 


i=i 




E hf E, WEf hf > 


i=i 


'ibk 

Skii-Pk) 


(19) 


'fc- 


TABLE I 

AO-based non-robust transceiver design algorithm 


1. Initialize W and define the tolerance of accuracy S. 

2. Repeat 

2.1 Solve problem GD with fixed W to obtain the updated 

{ffc 5 Pk }- 

2.2 Solve problem l |2Q) with fixed {f/-, pk} to obtain the up¬ 
dated W. Employ the proposed rank-one recovery method 
in Appendix if higher-rank solutions are returned by 
solving problem i |20^ . 

3. Until the total power consumption between two adjacent itera¬ 
tions is less than S or the total power consumption of step 2.2 
is higher than that of step 2.1. 


Hence, ( [T3] l can be reformulated as the following SDP problem 
by employing the well-known SDR technique 


K 

min X; 


W fc=l 

s.t. (18) and (|19 


+ Tr(GfcWG 
W A 0, V/fc G /C, 

3 Ijv • Different from 


Tr(W) 


( 20 ) 


where = (Gf^)^ 0 Ijv,.- Different from ( |15| l, the opti¬ 
mal solution W* to problem (20i is not necessarily rank- 
one; hence a simple rank-one recovery method based on a 
randomization procedure is proposed in Appendix[A|to address 
this issue. 

The AO-based iterative algorithm to solve problem ( [T3] ) is 
summarized in Table U Note that in each iteration of this 
algorithm, the objective function can only be decreasecf] and 
it is also lower bounded by zero, thus the convergence of the 
iterative algorithm is guaranteed. 


B. Robust AO-Based Joint Transceiver Design 

In this subsection, we address the robust counterpart of the 
non-robust AO-based joint transceiver design by employing 
the channel error model considered in Subsection III-BI We 
demonstrate that the corresponding subproblems can be con¬ 
verted to alternative forms where the concepts of SDR and 
S-procedure can be applied. 

First, let us consider the optimization of {{k,Pk} in (141 
when W is fixed under imperfect CSI , i.e. -f 

Gfc, V/c G JC. In this case, the SINR constraints in (15 i can be 
replaced by the following quadratic forms: 


OJt 


(hfe -f Bk) Ufc(hfc -f Bk) > -|-, 

Pk 

where Bk satisfy ED and Ufe can be expressed a; 


Ur 


1 ^ 

= Qfc y Qj 


cr'WW 


H 


j¥=k 


( 21 ) 


( 22 ) 


Similarly, the EH constraints can be reformulated as the 
following expression: 


(hfc+efe) Vfe(hfc-f efc)-f cTfc > 


i/’fc 


^/c(l - Pk) 


M\^kr<pi 

(23) 


^Due to the rank-one recovery method, the total transmission power of step 
2.2 is not necessarily less than that of step 2.1. However, we can terminate the 
algorithm if that happens to ensure the convergence of the algorithm. Note 
that this rare situation does not happen in all our simulations. 
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K 


where ^ 

j=i 


(T^WW 


H 


By applying the S-procedure, the constraints in (21 1 and 


can be reformulated as finite convex constraints, which 
are equivalent to the following two linear matrix inequality 
(LMl) constraints 


U/c + Afcl 

hfU, 


Ufch 




- ujlpk - Xkvl 


y 0 , 


(24) 


Vfc + /Tfel 


hf Vfehfe + cr^ - ^Qk - fJ-kvl 


h 0 , 


(25) 

where pk = qk = tz —^ while Afc > 0 and /i^ > 0 are 

Pk J- Pk 

slack variables. 


With (24 1 and (25 1 , the robust counterpart of problem ( [T5| ) 
can be reformulated as follows by ignoring the rank-one 
constraints for all F^^s: 

min X; Tr(Ffc)-f X; Tr(Qfe)-f cr^llWf 

{Ffc, pfc, gfc, Afc, /Xfc } • 

s.t. (24| and (251, Afe > 0, pk > 0, 

Pk Qk invp(pfe) -f invp(gfe) < 1, 

Ffc ^ 0, Vfc G 1C. 

(26) 

It is worth noting that different from problem which 
always returns a rank-one solution, the optimal solution {F^} 
of ( [26] l may not be of rank-one, in which case, additional 
processing steps may be needed to extract a rank-one solution 
from the F^. To this end, a rank-one recovery method based 
on a randomization procedure but with lower complexity is 
presented in Appendix [A| 

Next, we consider the optimization of the RS AF matrix W 
with fixed {fk,Pk}- Similar to Subsection III-A we have the 
following expression 

gkj = vec(conj(hfc)fjG'^) = (Gf^- (g) l 7 vJconj(hfc) (27) 


K 


(28) 


Ufch 


fell/c 


U/c ^ 

hfUfe hfUfchfc-cr^-w^Pfc-Afep^ 


^ 0, 


TABLE II 

AO-based robust transceiver design algorithm 


1. Initialize W and define the tolerance of accuracy S. 

2. Repeat 

2.1 Solve problem j26) with fixed W to obtain the updated 
{Ffc, pfc}. If (Ffcl are rank-one matrices, use their prin¬ 
cipal eigenvectors as a solution. Otherwise, employ the 
rank-one recovery method to recover a rank-one solution 
from {Ffc}. 

2.2 Solve problem ID with fixed {Ffc, pfc} to obtain the 
updated W. If W is a rank-one matrix, use its principal 
eigenvector as a solution. Otherwise, employ the method 
in Appendix [a| to recover a rank-one solution from W. 

3. Until the total power consumption between two adjacent itera¬ 
tions is less than <5 or increased. 


_ _ K _ Nr ^ 

where ^ W.-a^ ^ E^WEf and A^ > 0. 

i/fc fc=i 

Similarly, the robust version of ([T9| can be formulated as 


Vfc -f /ifel 


hfv. 


hFVrh 




- '^qk - Pkvl 


>■ 0 , 


K 


Nr 


(30) 


where Vk = Wfc 
f=i 


S EfcWE^ and pk > 0. Hence, 

fc=i 

problem (14i with fixed {ik, Pk} can be formulated as the 
following SDP problem with the SDR technique 


- Tr(CfeWG 


_ min Y} (iffe 
W, {.Wfc} k=i ^ _ 

s.t. (29) and (Is^, W A 0, Afc > 0, pfe > 0, Vfc, 


Tr(W) 


(31) 

where Cfc has been defined in Subsection |III-A| Note that the 
optimal solution of (31 1 is also not necessarily rank-one and 


= Gjconj(hfc), 

where Gj = Gfj (g) j G 1C. Thus, the robust version of 
(fTSll can be formulated as 


i(hf+ ef)Wfc(hfc + efc) 


- E (hf + ef )Wj (hfe -f Gfc) 

6 (hf + ef )E,WEf (hfc + efc) 

i=i 

^ kG 1C, 


where Wj = G![^Wconj(Gj), and W has already been 
defined in Subsection |III-A| By applying the S-procedure, we 
can transform (28 1 into the following LMIs 


(29) 


thus we can employ the method in Appendix [A| to address this 
issue. 

We summarize the AO-based robust joint transceiver design 
algorithm in Table For practical implementation, the opti¬ 
mally designed AF matrix W and PS ratios {pk} should be 
fed forward from the BS to the RS and receivers, respectively, 
through signaling channels prior to data transmission. 

IV. Switched Relaying Based Joint Transceiver 
Design 

In the previous section, we proposed novel AO-based 
transceiver design algorithms. As will be shown in Section 
VI these algorithms can achieve an excellent performance in 
transmission, but they are characterized by higher signaling 
overhead and computational complexity. In this section, mo¬ 
tivated by these considerations, we develop alternative SR- 
based transceiver design algorithms that are more efficient 
and simpler to implement. As illustrated in Fig. we equip 
the RS with a finite codebook of permutation matricej^ i.e., 
T = {Ti,T 2 ,...,Tb}, where T; G {0, is the l-th 

matrix in the codebook, index I G {1,... ,B} and B denotes 

^A permutation matrix is a square binary matrix that has exactly one entry 
equal to 1 in each row and each column, while all the other entries are equal 
to 0. 
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the codebook size which satishes B In order to 

reduce the signaling overhead and the number of optimization 
variables, the RS AF matrix W is constructed by multiplying 
the appropriate permutation matrix from the codebook with 
a power scaling factor. That is to say, the optimization of 
W is replaced by y^Ti, where Pi is a variable power 
scaling factor. Thus, each permutation matrix gives rise to 
a permuted channel matrix for which we can design a so- 
called latent transceiver, consisting of the BS beamforming 
vectors, the RS power scaling factor and the receiver PS 
ratios. Among the B latent transceivers so designed, by an 
algorithm to be developed below, the optimal one with index 
lopt is chosen by means of a suitable selection criterion for 
transmission. Specihcally, the selection mechanism is designed 
to choose the optimal latent transceiver with the minimum 
power consumption, which can be expressed as 


lopt = argmin Pi, 


(32) 


where Pi denotes the total transmission power corresponding 
to the Zth latent transceiver. 

Before data transmission, the BS only sends the index of 
the permutation matrix, the RS power scaling factor and the 
receiver PS ratios corresponding to the optimum transceiver 
to the RS and the receivers. Compared to the AO-based 
algorithms, the SR-based algorithms can signihcantly reduce 
the number of signaling bits and the computational complexity. 
The proposed SR-based scheme works as follows. 


• The BS designs the B latent transceivers based on avail¬ 
able permutation matrices within the codebook; it then 
determines the optimal latent transceiver based on (32i. 

• The BS sends the optimal permutation index lopt, the 
corresponding RS power scaling factor, and the PS ratios 
corresponding to the optimal transceiver to the RS and 
the receivers through signaling channels. 

• The RS constructs the optimal RS AF matrix based on 
the forwarded RS power scaling factor and the stored 
permutation matrix with index lopt- 

We introduce the non-robust and robust latent transceiver 
design algorithms in Subsection IV-A and IV-B[ respectively. 
In Subsection IV-C[ a simplihed SR-based transceiver design 
algorithm is proposed. The design approach for the codebook 
of permutation matrices is presented in Subsection |IV-D| 


A. Non-Robust Latent Transceiver Design 

Firstly, we introduce a non-robust algorithm to construct 
the kh latent transceiver. We aim to design the beamforming 
vectors {f^}, the RS power scaling factor Pi and the receiver 
PS ratios {p\} so as to minimize the sum of BS and RS 
transmit power under both SINR and EH constraints. We note 
that superscripts I in {f^}, {p^}, etc., have been removed in the 

®The total number of permutation matrices at the RS is Nr\. It is not 
realistic to use all the permutation matrices as the codebook when Nr is large. 
Thus, we propose to construct a codebook of permutation matrices with B 
elements, B <C NA-. The codebook design approach will be introduced in 
Subsection |IV-D| In practice, the value of B should be chosen to achieve 
a suitable tradeoff between performance requirements and implementation 
complexity. 


following for notational simplicity. For each latent transceiver, 
the optimization problem can be formulated as 


min Pg 

{ffc. Pk}, Pi 


Pi 


R 


s.t. Ffc > 7 fe, P^^ > t/ife, 
0<pfe<l, W = v^T,, Vfce/C. 


(33) 


As we can see, although we replaced the AF matrix W with 
^/P~lTl, problem (331 is still non-convex and difficult to solve 
due to the coupling between variables f^, pk and Pi. Our 
proposed method is motivated by the observation that problem 
(|^ can be reformulated as a difference of convex (DC) 
programming problem with proper transformations. Thus, the 
concept of CCCP |251, can be adopted to iteratively solve 


the DC problem, as explained below. 

First, we introduce the variable substitution 

1 


‘Pi = 


/3r 


(34) 


and further dehne the following vectors: 

p = [pi,---,pk]’^, q= [qi,---,qKr, 

where pk and have already been introduced in Subsection 


(35) 


III-B alongside with ([24)i and (25 1 . Then, the objective func 


tion of problem (33 i can be transformed into 


K 




(36) 


which is strictly jointly convex in the variables f} G IR+ x 

(|22j. 

Similar to the transformation of the objective function, the 
SINR constraints can be rewritten as 


Wfe(r) - Xk[r) < 0, 

where Wfc(r) and Xk{v) are dehned as 


(37) 


(38) 


K 

Wk{r) = ^ff G^TfhfchfTzGfj 

jpk 

+ hfc -f ^ujlipk + Plf, 

Xkir) = ]ujUpk - pif + -ff G^TfhfchfT,Gffc. (39) 

4 7fe 

Moreover, the EH constraints can be recast as 

Vkij) - Zk(r) <Q, (40) 

where 

yfc(r) = + Pif - hfc, (41) 


K 


Zfc(r) = ^ ff G^Tf hfehf TiGf, (42) 

i=i 

+ - Pif- 

We remark that (1^, (|g and are all convex 

functions jointly with respect to the variables in r G x 




















Beamforming 



Fig. 2. 


Proposed transceiver scheme with switched relaying processing for multiuser MISO relay system. The dashed line denotes signaling channels. 




X IR+ X Thus, problem (33 i can be equivalently 


reformulated as the following DC program: 


min P;(r) 

s.t. Wkir) - Xkir) < 0, yk{v) - zuir) < 0, (43) 

Pk > 1, <7fc > 1, invp(pfc) +invp(gfc) < 1, V/c S K, 


where the last set of inequality constraints must be satisfied 
with equality at optimality, for otherwise, the objective value 
can be further decreased by decreasing the pfc’s. 

According to the concept of CCCP, we approximate the 
functions aife(r) and Zfe(r) in the fth iteration by their first- 
order Taylor expansions around the current point denoted 
as and Zfe(r(®\r), respectively. With the help of 

and 1^, is given by 

Xfc(rW,r) =Xfe(rW)-f 27^{V:rfc(rW)^(r-r«)}, (44) 


where Va;fe(r^*)) denotes the conjugate derivative of the func¬ 
tion Xk (r) with respect to the complex vector rQWe note that 
irfc(r^*\r) is an affine function of r. is given by 


Va;fc(rW) = [Oix(fc-i), 0ix(2tf-fc), 




(45) 


Similarly, can be expressed as 

Zfc(r«,r) =Zfc(r«) + 27^{Vzfc(r«)^(r-r«)}, (46) 


where 


V^fe(rW) = [Oix(if+fe-i), ^(4*^ - Oix(K-k), 


irr2 _ tr'. tr', 

(9/c ‘Pi j) ) , . . . , (IjfcT^ j J , 

Gk = G^Tf hfehf TzG. 

(47) 

Then, in the ith iteration of the proposed CCCP based 
algorithm, we have the following convex optimization problem 


tk(Si) 


r{i)\T 


F(*)^TlT 


min P;(r) 

s.t. Wk[r) - Xk{r^"\r) < 0, yfc(r) - Zfe(r(*\r) < 0, 

Pk > 1, gfe > 1, invp(pfc) -f invp(( 7 fc) < 1, Vfc e /C, 

(48) 

whose solution is denoted by 


^ Since r is composed of both re al an d complex variables, we make a 
modification of such that j44| holds for r. 


TABLE III 

Non-robust latent transceiver design algorithm 


1. Define the tolerance of accuracy 5 and the maximum number 
of iteration Nmax- 

2. For I = 1,..., B 

3. - Initialize the algorithm with a feasible point which is 

obtained by solving problem \33) with W = T^. Set the 
iteration number i = 0. 

4. - Repeat 

- Compute the affine approxima tion and 

Sfc(r^^\r) acc ordin g to |44| and l |46) , respectively. 

- Solve problem j49t , and assign the solution to 

- Update the iteration number : i = i 1. 

- Until |Pi(r(^+i)) - Pi{r(^'>)\ < 5 or the maximum 
number of iterations is reached, i.e., i > A^max- 

5. - Obtain the Ah latent transceiver {f^, 

6. End 


Proposition 1: By introducing a new set of variables dk. 


dk, Ck and e^, k G JC, problem (48i can be reformulated as 
the following SOCP problem 


s.t. 



Pi+ P 2 + P 3 

.,f|,(Pi-l)/2]||<(Pi + l)/2, 

{ GiKr ,{ P2 - Pi ) m \\<{ P2 + pi )/ 2 , 


Pk > 1, gfe > 1, invp(pfc) -f invp(gfe) < 1, VA: e 1C, 


(49) 


where Wk is defined in ( |7T] l. 

Proof: Please refer to Appendix 

We can show that the proposed CCCP based iterative 
algorithm for the hh transceiver design converges to a local 
optimal solution of problem ( [33] ). The proof is similar to 
that of Lemma 2 and Theorem 1 in and we therefore 
omit the details. We summarize the proposed non-robust latent 


transceiver design algorithm in Table III 


B. Robust Latent Transceiver Design 

Secondly, let us propose a robust latent transceiver design 
algorithm by taking the channel errors into consideration. As 
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we can see, the concept of SR can also be applied to problem 
•HI, which can be reformulated as follows: 


min Pb 

{ffc, Pk}, Pi 

s.t. Tk > 7fc, 
W = yMTl 


Pr 


pEH 


>ipk, 0 < Pk < 1, 
efelP < vl, V/c e K.. 


(50) 


We employ a subgradient-type iterative algorithm to solve this 
problem. Our proposed method is motivated by the observation 
that, if we fix /?;, then problem ( [50| is equivalent to problem 
(26 1 with the SDR techniqueWith fixed Pi, (24i and (251 
can be reformulated as the following two LMIs: 


Ufc -t- Afcl 




Vfchfe 


^ 0 , 


Vfc -f ^fcl 

h»V, 


Sf v.hk + If - ^ 

K ^ K 


(51) 


(52) 


where Qj - k = Qj +^^1, and 






Qfc = TiGFkG^Tf. 


Then, problem (|5^ with the use of the SDR technique and 
fixed Pi can be expressed as 


K 


m) = 

min X) Tr(Ffc) + pi Tr(Qfc) -f a'^Pi 

{Ffc, Pfc, Qki Mfcj k—1 k—1 

s.t. (51) and (52), Xk > 0, > 0, 


Pk > 1, > 1, invp(pfc) + invp(gfc) < 1, 

Fk h 0, Vfc e /C. 

(53) 

We here consider the SDR version of problem ( [50l l instead of 
its original form with rank-one constraints, since strong duality 
holds if problem ( [5^ is feasible. The partial dual problem of 
([5^ can be formulated as (5^ shown at the top of next page. 


where Xj. = 


Xfe 


Xfc 


Yfc 

yf 


Yfc 


Ik 

15T1) and 


denote 


and Yfc = 

_ Xk j _ 

the dual variables associated with constraints 
respectively. It is worth noting that X^ and Y k can be obtained 
as side information provided with any standard SDP solver, 
since dual variables are served as a certificate for optimality. 

Next, we can use the subgradient method pO) to iteratively 
solve the SDR version of problem (50 1 . At the (i -f l)th 
iteration, the power scaling factor Pi can be updated according 
to 


Pi{i + l) = [Pi{i)-9{i)s{i)] + , (55) 

where [•]f = max (.,£), 9{i) is the step-size in the ith 
iteration and s{i) denotes a subgradient of /(/?/) at Pi{i). The 


^The concept of CCCP can also be applied to the robust problem iSOj. 
However, based on simulations we find that the performance of the CCCP 
based algorithm is inferior to that of the subgradient-type algorithm. Thus, 
we employ the subgradient-type algorithm instead. 


TABLE IV 

Robust latent transceiver design algorithm 


1. Initialize the step-size 9, define the tolerance of accuracy 5 and 
the maximum number of iterations Nmax- 

2. For / = 1,..., B 

3. - Initialize ft(0). Set the iteration number 2 = 0. 

4. - Repeat 

- Solve problem j53| to obtain the optimal and 

the con'esponding dual variables 

- Update the power scaling factor according to (55). 

- Update the iteration number : 2 = 2 -|- 1. 

- Until |/(ft (2 -|-1)) —/(ft( 2 ))| < 5 or the maximum 
number of iterations is reached, i.e., 2 > A^max- 

5. - Obtain the hh latent transceiver {f^, pj, 

6. End 


subgradient s{i) can be calculated as ||^ 


K if 2 2 

cri UJ^ 


s(*) = 


+ (56) 


fc=l 

K 

^Vki 

k^l 


k=i pWii^ 

'P’k 




where = T/GF^G^T^, {F^,pU denotes the optimal 
solution of problem ( [5^ and {x]., yl} denote the lower right 
comer elements of X^ and Y^, which are the optimal dual 
variables associated with (|5T| and ( |52| ). 

Finally, the robust late nt transceiver design algorithm is 
summarized in Table m 


C. Proposed Simplified SR-based Transceiver Design 

the proposed 


As shown in Subsection IV-A and IV-B 


SR-based transceiver design algorithm involves devising B 
latent transceivers corresponding to the elements in T where 
for each transceiver, we employ iterative methods to address 
the highly non-convex problem. Specifically, the design of 
each transceiver involves an iterative algorithm to obtain the 
corresponding BS beamforming vectors, RS power scaling 
factor and receiver PS ratios. However, for given CSIs only the 
best transceiver with the minimum total power consumption is 
selected for transmission according to the selection mechanism 
( |32| ) while the other B—\ solutions are discarded. Thus, there 
will be a considerable waste in the computational resources. 
In order to improve the proposed SR-based transceiver design 
algorithm and make it more suitable for practical implemen¬ 
tation, we propose a heuristic approach, referred to as the 
simplified latent transceiver design algorithm, to address the 
aforementioned problem. Based on simulation experiments, 
we find that a good initial point obtained by solving problem 


or (501 with fixed initial power scaling factor Pi almost 
always leads to a better convergent point than by solving 


®For the non-robust case, we can also employ the subgradient method to 
address problem with each subproblem being an SDP problem. It is well 
known that solving an SDP problem requires relatively high computational 
complexity compared with solving an SOCP problem. Moreover, it is impor¬ 
tant to note that the two algorithms for the non-robust case have very close 
performance in our simulations. Thus, we employ the CCCP based algorithm 
for the non-robust case. 
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K 

f(6i) = max min — Y] Tr 
{X,,Ya {Fk,Pk} /Yl 


Ufc + Afcl 
hfUfc 



k '-’k 

Vfc + //fel Vfchfc 

, hfV,hfc + f 


Afcjy^ 


hfV, 


Ufchfc 

Pi 

E Tr(Ffc) 

fc = l 


■Pi f: Tr(Qfc) +1x2/3,, 

A;=l 


(54) 


TABLE V 

Simplified non-robust latent transceiver design algorithm 


smallest singular value, i.e., 


1. Define the tolerance of accuracy (5 and the maximum iteration 
number Nmax- 

2. For the hh latent transceiver (/ = 1,..., S) 

3. - Compute a initial feasible point which is obtained by 

solving problem j33t with W = T^. 

4. end 

5. Choose the initial and ^ with the minimum objective 
value. Set the iteration number i = 0. 

6. Repeat 

- Compute the affine approximation and 

r) acc ordi ng to |44| and l |46| , respectively. 

- Solve problem j49t , and assign the solution to 

- Update the iteration number : i = i 1. 

1. Until |Pi„^*{r(»+i)) - P,„^,(r«)| < 5 or the maximum 
number of iterations is reached, i.e., i > A^max- 


with a not so good initial point. Thus, we only design the 
particular transceiver with the best initial point. The proposed 
simplified non-robust latent transceiver design algorithm is 
summarized in Table [V] The robust version of the simplified 
latent transceiver design algorithm is similar to the non-robust 
case, and thus is omitted here. We will show in Section [Vl] 
that the simplified SR-based transceiver design algorithms can 
achieve a similar performance with that of the proposed SR- 
based algorithms in Table |III] and [rV| 


D. Codebook Design 

In this subsection, we propose two algorithms to construct 
the codebook of permutation matrices. The basic principle of 
the proposed algorithms is to choose the permutation matrices 
which are more likely to result in lower transmission power. 
We make a promising observation that when the singular 
values of the permuted channel matrix HT/G are larger 
or more equally distributed, the total power consumption is 
usually smaller, where the channel matrix H has already 


been defined in Subsection II-C and HT,G represents the 
equivalent channel matrix between the BS and the receivers. 

Let 7r[ denotes the A:th singular value of HT,G, k G 1C. 
The first scheme (referred to as Sum-Max method) constructs 
the codebook of permutation matrices by choosing the ones 
which correspond to the B largest sums of singular values. 


I.e., 


{Ti,...,Tb} = argmaxB 

Ti 



(57) 


where maxB( ) returns the permutation matrices correspond 
to the B largest values of its argument. 

Alternatively, the second scheme (referred to as Max-Min 
method) chooses the permutation matrices to maximize the 


{Ti,..., Tb} = argmaxB [ min(7r[) j . (58) 

It is worth noting that the proposed two schemes are heuristic 
techniques which may serve as a shortcut to the process 
of finding a satisfactory solution. The performance of the 
codebook design algorithms will be studied in the simulation 
results in Section Iyj 


V. Complexity Analysis 


In Section m and HYl we proposed the AO-based and 
SR-based transceiver design algorithms for problem ( [T3| ) and 
respectively. In this section, we compare the relative 
computational complexity of the proposed transceiver design 
algorithms. Moreover, we apply the same basic element of 
complexity analysis as in p2) . Among the proposed algo¬ 
rithms in this paper, we consider: 

• Non-robust AO: the algorithm summarized in Table |I] 

• Robust AO: the algorithm summarized in Table [I^ 

• Non-robust SR: the algorithm summarized in Table 

• Robust SR: the algorithm summarized in Table 

• Simplified non-robust SR: the algorithm summarized in 
Table 0 

• Simplified robust SR: the robust counterpart of the algo¬ 
rithm in Table |V] 


The complexity of the non-robust AO algorithm is domi¬ 
nated by solving problems p5] ) and ( p0| ) Ii times , where Ii 
denotes the iteration number. We note that the dual problems 
of ( [TSl l and ( |20l i can be solved instead of ( [T5| ) and ( |20| ) for 
better efficiency. Consider the dual problem of (15i, which 
involves K LMI constraints of size Nt and on the order of 
ni = 0{KN^ -\- K) decision variables. Thus, the complexity 
of a generic interior-point method for solving problem 0 is 
given by 0(ni\/KNt{KNf riiKNf -|- nf)). Similarly, the 
complexity of solving the dual problem of ( [20l l can be written 
as + n 2 N^ -\- n|)), where n 2 = 0{NE is the 

order of the corresponding decision variables. The complexity 
of solving problem ( [65] l can be neglected since it admits 
closed-form solution. Thus, the overall complexity of the non- 
robust AO algorithm is on the order of the quantity shown in 
the first row of Table m 

The complexity of the robust AO algorithm is dominated by 
solving problems (261 and ( [3T] ) I 2 times, plus the complex¬ 
ity of the rank-one recovery method. Problem ([26| involves 


KNf AK variables, 2K LMI constraints of size Nr A- I 
and K LMI constraints of size Nt- Problem iB involves 
Nf -I- 2K variables, 2K LMI constraints of size Nr 1 and 
1 LMI constraint of size Nf. Moreover, the complexity of 
the rank-one recovery method in Table VII is dominated by 
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solving problem ( |64l) R times, where R is the number of 
randomization steps|‘^| Problem (64 1 involves 2K variables 
and 2K second-order cone (SOC) constraints of dimension 
3. Thus, the overall complexity of the robust AO algorithm is 
on the order of the quantity shown in the second row of Table 

Iyj 

The complexity of the non-robust SR algorithm is dominated 
by solving problem ( |49| ) BI 3 times, where I 3 is the iteration 
number, since the complexity of computing the affine approx¬ 
imation and Zk{r^''\r) is negligible compared to 

solving ( |49] l. Problem ( |49] l involves 2K -f 3 SOC constraints, 
including 1 SOC of dimension KNt + 2, 1 SOC of dimension 
KNj. + 2, K SOCs of dimension K + 2, and K + 1 SOCs 
of dimension 3. The number of variables is on the order of 
0{KNt -f 2K). It follows that the complexity of the non- 
robust SR algorithm is on the order of the quantity shown in 
the third row of Table 

The complexity of the robust SR, simplified non-robust SR 
and simplified robust SR algorithms can be analyzed in a 
similar way; the corresponding complexity figures are shown 
in the fourth to sixth rows of Table |VI| respectively, where 
I 4 denotes the iteration number of the (simplified) robust 
SR algorithm. It is of interest to investigate the asymptotic 
complexity of the proposed algorithms when Nt, Nr and K 
are large, i.e., when we let Nr = Nt = K ^ 00 . We 
further assume that Ii = I 2 = I 3 = I 4 = I for simplicity. 
Under these conditions, one can verify that the complexities 
of the proposed algorithms in Table |Vl] are on the orders of 
21N^^, 2^1 Nl^, QBINl^, 3y/3BIN^°, 6 {B + I)Nf-^ and 
3\/3{B -f I)Nt^, respectively. As seen, the robust algorithms 
always consume more computational resources than their 
non-robust counterparts, the SR-based algorithms have lower 
complexity compared with the AO-based algorithms and the 
simplified SR-based algorithms have the lowest complexity. 


VI. Simulation Results 

In order to evaluate the performance of the proposed 
transceiver designs, numerical results have been obtained 
by performing computer simulations. In the simulations, we 
assume that both the first phase and the second phase channel 
coefficients are flat-fading i.i.d. with unit variance Rayleigh 
distribution. The nominal system configuration is defined by 
the following choice of parameters: Nt = Nr = K = 3, 
^ = 1, = —30dBm, uij. = uP' = —20dBm, 

and (T^ = — 30dBm unless otherwise specified. In addition, 
we assume equal SINR and EH thresholds at the destination 
receivers, i.e., 7^ = 7, 4’k = ip^ ^ ^nd equal 
norm bounds for the second phase channel error vectors, i.e., 
Vkj = V, Vj, k for simplicity. In the implementations of 
the various algorithms, the tolerance parameter is chosen as 
6 = 2x 10-3 while R = R = 20 and R = R = 50; all the 
rank-one recovery methods employ R = 100 randomization 
steps. All convex problems are solved by CVX mi on a 
desktop Intel (13-2100) CPU running at 3.1GHz with 4GB 
RAM. 


In the AO-based transceiver design algorithms, the RS AF 
matrices W can be initialized in three different ways: 

• Init-1 : Initialize the RS AF matrix as an identity matrix. 

• Init-2 : Random matrices are generated according to 
normal distribution with zero mean and unit variance. 

• Init-3 : According to Section III and IV the SR-based 
transceiver can be used to initialize the AO-based designs. 

Fig. shows the average power consumption performance 
versus the receiver EH target ip for the AO-based transceiver 
designs with the three different initialization methods. From 
the results, we can see that the best performance is achieved 
by Init-3 for both the non-robust and robust algorithms; Init-1 
and Init-2 result into similar performance. For example, the 
power consumption of Init-3 is 0.5dB less than that of Init- 
1 and Init-2 for the non-robust case, and the power saving 
rises up to IdB for the robust case. As seen, the AO-based 
algorithms are sensitive to the initial point. The use of Init-3 
is not realistic in practice since it has to employ both the AO- 
based and the SR-based algorithms. However, in the following 
simulations, we employ Init-3 as the initialization method for 
both the non-robust and robust AO algorithms, which serves 
as a performance bound unless otherwise stated. 

Fig. 1^ shows the performance comparison of the SR-based 


transceiver design algorithms proposed in Subsection IV-A 


and |IV-B| and the simplified SR-based transceiver design 
proposed in Subsection |IV-C| It is worth noting that the 
algorithms with codebook size B = 1 only use the identity 
matrix as the permutation matrix (naive method|^ and the 
algorithms with the B = 8 codebook employ the Sum-Max 
codebook design method. We include the identity matrix into 
the B = 8 codebook regardless of the singular values in 
order to guarantee that the performance of the algorithms with 
B = 8 is always better than with B = 1. We can see that 
the simplified SR-based transceiver design algorithms achieve 
almost the same performance as the non-simplified ones for 
both non-robust and robust cases. However, the simplified SR- 
based algorithms consume much less computational resources 
compared with their non-simplified counterparts. Thus, we 
employ the simplified SR-based design in the following sim¬ 
ulations for comparison unless otherwise specified. 

We compare the performance of the codebook design ap¬ 
proaches, i.e., the Sum-Max and Max-Min methods, in terms 
of average power consumption. For completeness, we also 
consider the performance of a codebook whose elements (per¬ 
mutation matrices) are randomly selected. The results which 
are illustrated in Fig. show the average power consumption 
performance curves versus EH target. We can see that the Sum- 
Max method outperforms other methods. Compared with the 
randomly generated codebook, the proposed Sum-Max method 
can lead to a power saving of 0.4dB for the non-robust case 
and 2dB for the robust case. The performance of the Max-Min 
method is slightly inferior to the Sum-Max method. 

Fig. 0 shows the average power consumption performance 
for various codebook sizes in the SR-based transceiver design 
algorithms. It is observed that the SR-based robust transceiver 


^^The complexity to recover W* is negligible since problem l |65| has 
closed-form solution. 


^^The naive method means that the RS only processes the received data 
vector Yfii by adjusting its power (i.e. scaling) and forwarding it. 
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TABLE VI 

Complexity analysis of the proposed transceiver designs 


Proposed design 

Complexity Order (suppressing the ln(l/£:) ) 

Non-robust AO 

hiOim^KNtiKN^^ + niKNf + nf)) + 0{n2Nr(N^ + n 2 N^ + n|))), m = OiKNf + K), na = 0(N^} 

Robust AO 

l2{0(n^^2K{Nr + 1) + KNt(2K(Nr + 1)“ -b KNf + 2n2K{Nr + 1)^ + n^KN'^ + nl))+ 
0{n2yj2K(Nr + \) + N^(2K{Nr + 1)3 + N^ + 2n2K{Nr + 1)^ -1 -b nl))+ 

R(0{2KVIK(2K?P' + 4^-2)))), m = 0(KNf + 4K), na = 0{N* + 2K) 

Non-robust SR 

Bl3(0{nV4K + 6i{KNt + 2)2 -b (KNr + 2)2 -b K{K + 2f + {K + 1)32 -b n^))), n = 0(KNt + 2K) 

Robust SR 

Bl4,{0{n3y/2K(Nr + 1) -b KNt(2K{Nr + 1)^ + KNf + 2niK(Nr + if + n^RNf + nf))) 

R{0{2KVIK{2K3'^ + 4X2)))), ^ 0(KN^ + 4K). 

Simplified non-robust SR 

{B + l3){0(nV4K + e({KNt + 2)2 -b {KNr + 2)2 + K{K + 2)2 + (X -b 1)32 + n2))), n = 0{KNt + 2K) 

Simplified robust SR 

(B + l4,)(0(n3^2K(Nr + 1) + KNt(2K{Nr + 1)^ -b KNf + 2niK(Nr + if + niKNf + nj))) 

-hX(C>( 2XV4X(2X32 -b 4X2)))), ^ 0{KN^ + 4X). 



Fig. 3. Initialization methods comparison of the AO-based transceiver designs 
(?? = 0.1, 7 = lOdB). 



Fig. 4. Performance comparison of the non-simplified SR-based transceiver 
designs and the simplified SR-based transceiver designs (77 = 0.1, 7 = lOdB, 

A(0) = 1,6 = 6). 


design algorithm with codebooks of size B = 4, 8,16 and 24 
achieves a power saving of 5dB compared to that with the B = 
1 codebook, and its non-robust counterpart with codebooks of 
size i? = 4,8,16 and 24 achieves about IdB in power saving 
compared with the B — 1 codebook|^We can also see from 

*^Note that the algorithms with S = 24 codebook conduct a search for 
all the Nr\ permutation matrices, which is known as the optimal SR-based 
transceiver design algorithm. 



Fig. 5. Comparison of the codebook design methods for the SR-based 
transceiver designs (77 = 0.1, 7 = lOdB, /3;(0) = 1, 6 = 6 ). 

this figure that an increase in B does not necessarily gives 
rise to an improvement of the power consumption performance 
when B > 8, which is the reason we proposed to construct a 
restricted number of transceivers, i.e., to select a small fraction 
of all the permutation matrices for transceiver design (B 
Nrl). 



\|/ {dBm) 


Fig. 6 . Codebook size comparison of the SR-based transceiver designs (77 = 
0 . 1 , 7 = lOdB, ft( 0 ) = l,e = 6 ). 

In the next series of simulations, we examine the compar¬ 
ative performance of the AO-based and SR-based transceiver 
design algorithms. Fig. shows the average normalized trans- 
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mitted power versus ip for these two algorithms. It is observed 
that the performance of the SR-based algorithms is very close 
to that of the AO-based algorithms. From the complexity 
analysis in Section |VI[ we note that if increases, the 
complexity of solving problem (20i and •HD in the AO- 
based transceiver design might become unacceptable, which 
limits the practicality of the algorithm. Thus, the SR-based 
transceiver design is very promising and suitable for systems 
with large Nr- It is also important to mention that the number 
of signaling bits of the SR-based algorithms is much less 
than the AO-based algorithms since only the index of the 
optimal transceiver and the power scaling factor have to be 
sent to the RS instead of the whole RS AF matrix. Besides, 
a comparison of the feasibility rate between the AO-based 
and SR-based algorithms is provided. Specihcally, Fig. 
presents the feasibility rate comparison versus SINR target 
7. One can observe that the feasibility rate performance of 
the AO-based algorithms is slightly inferior to that of the SR- 
based algorithms since we employ Init-S as the initialization 
method for the AO-based algorithms. Fig. [^presents a similar 
comparison of the feasibility rate versus the channel error 
bound rj. Similarly, from this hgure, we can see that the 
performance of the AO-based algorithms and the SR-based 
algorithms is very close. The non-robust algorithms fail to 
satisfy both the SINR and EH constraints almost all the time 
under the NBE model. 



Fig. 7. Performance comparison of the SR-based transceiver designs and the 
AO-based transceiver designs (77 = 0.1, 7 = lOdB, A(0) = 1, 0 = 6 ). 



Y(dB) 


Fig. 8 . Feasibility rate comparison of the SR-based transceiver designs and 
the AO-based transceiver designs versus SINR threshold 7 (77 = 0.1, A (0) = 

l,e = 6). 
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Fig. 9. Feasibility rate comparison of the SR-based transceiver designs and 
the AO-based transceiver designs versus CSI error bound rj {-y = lOdB, 
ijj = 5dBm, /3j(0) = 1, 9 = 6 ). 


matrix codebook in the SR-based algorithms were proposed. 
The simulation results validated the effectiveness of the pro¬ 
posed joint transceiver design algorithms and the necessity of 
using a robust formalism in the presence of imperfect CSI. We 
note that the proposed algorithms can be easily extended to 
more complicated (multi-hop) relay systems, which remains 
an open avenue for future work. 


VII. Conclusion 

In this paper, we considered the joint transceiver design 
problem for multiuser MISO relay systems with energy har¬ 
vesting. We proposed AO-based and SR-based algorithms, 
including both non-robust and robust versions to CSI errors, 
for the joint optimization of the BS beamforming vectors, 
the RS AE matrix and the receiver PS ratios. Especially, 
the SR-based algorithms can achieve almost the same per¬ 
formance compared with the AO-based algorithms but with 
much reduced computational complexity and overhead. We 
also presented a simplihed SR-based algorithm and carried out 
a detailed complexity analysis of all the proposed algorithms. 
Two efficient approaches for the design of the permutation 


Appendix A 

Proposed Rank-one Recovery Method 

As we mentioned in Section [III and|^ the optimal solution 
of problem @, (|^ and (|53|l is not guaranteed to be 
rank-one. In this appendix, we propose a rank-one recovery 
method based on a randomization procedure. Inspired by the 
concept of worst-case robustness p^-ph), we consider the 
following subproblems: 

min |(hf+ ef)W*Gf*7 (59) 
\\ek\r<vt 

vuj = inax J(hf + ef )AV*Gf*7 j ^ k, (60) 
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||(hf-f ef)W*f, 

(61) 






(62) 

1. 

\\ ekP<vl 




where and W* are extracted from the optimal solutions 
and W* using the randomization procedure as will be 


detailed in Table 


VII 


It is not difficult to see that problem 
I can be solved by the Cauchy-Schwarz inequality and 
Lagrange multiplier method. We omit the detailed derivation 
and let denote the corresponding optimal 

objective values. 

For problem (26i and (53 i, we propose to scale up by 


and then jointly optimize and {pk} to recover 

a rank-one solution. The optimization problem can be formu¬ 
lated as 


TABLE VII 

Proposed rank-one recovery method 


For the ith randomization step {i = 1,..., R, R is the number 
of randomization steps.) 

- If i = 1, let equal to the principal component of 
F^. Else, calculate the eigenvalue decomposition = 
QAQ"^, then the candidate vector is generated as = 
QA2 w, where w is an i.i.d. complex Gaussian vector 
with zero mean and unit variance. 

- Solve problem i|64| to obtain the optimal P^.}- If ths 

problem turns out to be infeasible, discard that candidate 
vector. Otherwise, save the candidate {f^, and the 

objective value of the problem. 

End 

Select the candidate {f^,(pj,p^} which corresponds to the 
minimum objective value. 


K K 

min ^ 

{Vk,Pk} k=l k=l 

„ t _ PkVkulk 

5..L. ^ 

Pk H (pjvl.-\-pk(T^wl+pkcrl-\-(^l 
j¥^k 




> Ik, 


(63) 


K 


E + ^vK +(tI> 


i=i 




<Pfc > 0, 0 < pfc < 1, Vfc e /c. 

The above problem can be reformulated as the following 
SOCP problem 

K 


^k 'f:+ E ¥^fc||W*Gf^ 


min X; 

{‘Pfcr PkJ k —1 k — 1 

s.t. \\[ 2oJk,Zk - Pk]\\ < Zk+ Pk, 
ll[ 2 \/VE^:^fc - 1 + /Ofc]ll <Zk + l- Pk, 

K 


VkU^k 


Zk = - E ^j'^kj -^k- 

^ K ^ 

Zk= + ^rW*k + 




k'> 


min E ‘Fl|W*Gf*| 

{V,Pk} k=l 

^ __ PkvKk 


Appendix B 

The proof of Proposition 1 

We first introduce four sets of auxiliary variables dk, dk, Ck 
and ek,yk G 1C, which satisfy 

dk = cr^hf hfc -f (66) 

4 7fc 

dk = -alpi + f - ip["^)iPk - Pi) + —T^{ft'’^Gk{k}, 


K 


Ik 


efc --cr^n^. Lrfct^. + — (q^ - P; j , (68) 

4 = 1 


(67) 


K 


(64) efe = 


27?.{^ ^ - p|*^)(9fc - Pi) + 

4 = 1 




(69) 


Pfc > 0, 0 < pfc < 1, Vfc e /c. 

Similarly, we also propose to scale up W* by y4^ and then 
jointly optimize {^/p} and {pk} to recover a rank-one solution 
for problem (2^ and ( [3T] l. The optimization problem can be 
formulated af 3 ^ 


where dk and Ck are constant scalars while dk and Ck are 
affine functions in r. Then, the SINR constraints in problem 


(48 1 can be rewritten as the following SOCs 

rp dk dk 1 


'ic ) 


dk — dk + I 

< n 


(70) 


where 


(65) 


-a.Vl|W*„ 

-- > 7fc, 

PkV ’'^ki+VPk'^r'^l+Pk'^l+i^l 

p E + ^l> ottEt) ’ 

4 = 1 

P > 0, 0 < pfe < 1, Vfc G /c. 

As has been discussed in our previous work pT) , the above 
optimization problem admits a closed-form solution and we 
omit the detailed derivation in this work. The propo sed rank- 
one recovery method is summarized in Table VIl|*^ 


Wfc = [hfT,Gfi,...,hfT,Gffe_i,hfTzGffe+i, 
...,hfT,GfK,iwfc(pfe + P/)]^. ^ ^ 

Similarly, the EH constraints in problem ( |48] l can be reformu¬ 
lated as the following SOCs 


l/’fc / , N Cfc — Cfc — 1 

!&(“ + ■"')’- 2 




(72) 

Next, by introducing slack variables Pi, P 2 and P 3 , we can 


^For problem j53|, W* can be expressed as where p* denotes 

relay power scaling factor corresponding to the 1th latent transceiver. 


equivalently write the objective function P; (r) of problem (48 1 
as follows 

Pi(r)=Pi+P2+P3, (73) 


the relay power scaTmg factor corresponding i 

^^For problem j2o| , is the solution of j59|-|62| with 

G ^ — 0 ^ k . 

^^We only list the detailed steps to recover f^, which can be easily extended 
to the recovery of W*. 


where 


[fi^,...,f^,(Pi-l)/2]|l<(Pi + l)/2, 


(74) 


(Gfi)^,..., (Gf^)^, (P2 - P/)/2]|| < (P2 + Pi)/2, (75) 
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\\[2ar,ifi - Ps]]] < + Ps- (76) 

Therefore, problem ( [48] l can be shown to be equivalent to the 
following SOCP problem 


min Pi 
r, Pi , P2 , P3 


7^2 + Ta 


s.t. (1^, (1^, ^ and (1^, 

Pk > 1, g/c > 1, invp(pfe) + invp(gfe) < 1, VA: G AC. 

This completes the proof. 


(77) 
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